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Abstract 

We have calculated the one-body Fermi and Gamow-Teller charge-current, and vector and axial- 
vector neutral-current nuclear matrix elements in nucleon matter at densities of 0.08, 0.16 and 0.24 
fm“^ and proton fractions ranging from 0.2 to 0.5. The correlated states for nucleon matter are 
obtained by operating on Fermi-gas states by a symmetrized product of pair correlation operators 
determined from variational calculations with the Argonne vl8 and Urbana IX two- and three- 
nucleon interactions. The squares of the charge current matrix elements are found to be quenched 
by 20 to 25 % by the short-range correlations in nucleon matter. Most of the quenching is due 
to spin-isospin correlations induced by the pion exchange interactions which change the isospins 
and spins of the nucleons. A large part of it can be related to the probability for a spin up proton 
quasi-particle to be a bare spin up/down proton/neutron. Within the interval considered the 
charge current matrix elements do not have significant dependence on the matter density, proton 
fraction and magnitudes of nucleon momenta; however, they do depend upon momentum transfer. 
The neutral current matrix elements have a significant dependence on the proton fraction. We 
also calculate the matrix elements of the nuclear Hamiltonian in the same correlated basis. These 
provide relatively mild effective interactions which give the variational energies in the Hartree- 
Fock approximation. The calculated two-nucleon effective interaction describes the spin-isospin 
susceptibilities of nuclear and neutron matter fairly accurately. However > 3-body terms are 
necessary to reproduce the compressibility. Realistic calculations of weak interaction rates in 
nucleon matter can presumably be carried out using the effective operators and interactions studied 
here. All presented results use the simple 2-body cluster approximation to calculate the correlated 
basis matrix elements. This allows for a clear discussion of the physical effects in the effective 
operators and interactions. 

PACS numbers: 21.30.Fe, 23.40.Hc, 26.50.-hx 
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I. INTRODUCTION 


Weak interactions in nucleon matter occur during the beta-decay of nuclei, electron and 
muon capture reactions, neutrino-nucleus scattering and in various astrophysical environ¬ 
ments, such as evolving stars, neutron stars and supernovae. They have been studied since 
Fermi proposed the hrst theory of beta-decay in 1934. Recently there has been much inter¬ 
est in weak interactions in the sun § 0 . those of and due to their use in neutrino 
detectors searching for neutrino oscillations 10) B 0 ) and in interactions of neutrinos with 
dense matter in neutron stars and supernovae |]^. Low energy weak interactions proceed 
mainly via the nuclear matrix elements of the following four one-body operators: 
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Here i is the nucleon number label and q is the momentum given by the weak boson to the 
nucleon. The Fermi coupling constant multiplying these operators is omitted for brevity, 
dw is the electroweak mixing angle, and qa is the ratio of the weak axial vector and Fermi 
coupling constants of the nucleon. The four operators are called Fermi (F), Gamow-Teller 
(GT), neutral-vector (NV) and neutral-axial-vector (NA). In the nonrelativistic domain, 
neglecting weak pair currents, the interaction of low energy neutrinos with nuclei and nucleon 
matter and nuclear beta-decay rates are proportional to the square of the matrix elements 
of these operators between initial and hnal nuclear states. 

Due to the strong forces, nuclear wave functions are highly correlated p, |^, and it is 
difficult to calculate the nuclear matrix elements. Using quantum Monte Garlo and Faddeev 
methods to calculate nuclear wave functions from realistic models of nuclear forces, the 
beta-decay matrix elements have been calculated for light nuclei with A <7 |^, |^. The 
calculated values for ^H, ®Li and ^Be are within 5 % of the observed, and better agreement 
is obtained after including weak pair currents. The weak muon capture by ^He has also been 


calculated [12 with realistic wave functions with similar success. 
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However, complete many-body calculations are not yet possible for nuclei like and 
heavier, as well as for nucleon matter. Most studies of weak interactions in these systems 
use effective interactions and shell-model and Fermi-gas wave functions in hnite nuclei and 
nucleon matter respectively. The random phase approximation (RPA) is commonly used. 
The pioneering work on GT transitions has been reviewed by Arima et. al 0. Some of 
the recent works are: [|^, |K] in [^, in p/-shell, and [|^ in neutron stars and 
supernovae. Typically the calculated rate of weak interactions is larger than observed; for 
example, a factor of ~ 0.6 brings the calculated p/-shell GT transition rates in agreement 
with experiment. Recent LSND results of charged current reaction cross sections of z/g Q 
and 1^ on ^^G are lower than the theoretical expectations by up to 20 %. 

This is not surprising since effective operators which take into account the effects of short 
range correlations, and not the bare operators given by Eq. (1) to (4), must be used along 
with effective interactions as is well known from the works of Arima and collaborators [0 . 


In nuclei near the line of stability the observed spectra and beta-decay rates have been 
used to model the effective interactions and operators, but for neutron stars and supernovae 
matter we have to calculate them from realistic models of nuclear forces. In p/-shell and 


heavier nuclei the effective interaction is also obtained from bare forces 19 


There are several ways to obtain consistent sets of effective operators and interactions 
starting from a bare nuclear Hamiltonian. For example, one can introduce a model space 
and employ the Lee-Suzuki similarity transformation a as in the no core shell model type 
approach |0. In this theory the effective operators and interactions take into account 
the truncated Hilbert space. They are used in the retained model space to predict the 
observables. In the present work we use the correlated basis (GB) approach 0, 0, evolved 
out of variational theories of quantum liquids In this theory the uncorrelated shell model 
or Fermi-gas states are transformed by correlation operators to GB states without truncation 
of the Hilbert space. The effective operators and interactions are matrix elements of the bare 
quantities in the GB states; they take into account the effects of short range correlations. 
The correlation operators are chosen such that the nuclear interactions are relatively mild in 
the GB. Observables are calculated using standard many-body perturbation theory methods 
in GB. 

Here we focus on weak interactions in nucleon matter. In variational calculations IPI the 
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nuclear matter wave functions are approximated with correlated states: 




( 5 ) 


i<j 


where are uncorrelated Fermi-gas (FG) states and Fij are pair correlation operators. 
The iSn denotes a symmetrized product necessary because the Fij and Fik do not com¬ 
mute. One can also relate uncorrelated shell model states to correlated states in a similar 
way. The correlated states obtained from Eq. (|^) are not orthogonal; we assume that they 
are orthonormalized using a combination of Lowdin and Schmidt transformations pre¬ 
serving the diagonal matrix elements of the Hamiltonian. However, the orthonormalization 
corrections are of higher order than those considered here. 

Let |X) denote the orthonormal correlated states. The effective interactions in CB per¬ 
turbation theory are dehned such that: 


{X\H\Y) = {^x\Ho + Hi\^y) , 


h, = y. 
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Here H is the nuclear Hamiltonian containing realistic two- and possibly three-nucleon 
interactions. Even when if has only two-body interactions the CB FT/ can have three- and 
higher body terms. Since the correlated states are expected to be close to the eigenstates 
of H, the non-diagonal matrix elements {X ^ Y\FI\Y) are small. This implies that the 
CB effective interactions can be used in perturbation expansions based on the Hartree-Fock 
approximation. However, the 1st order results are often not sufficiently accurate. The 
product of pair correlation operators (Eq. (D) can not transform the uncorrelated states 
into the exact eigenstates of H. CB calculations of the optical potential of nucleons in 


nuclear matter including up to 2nd order terms in iF/, and of the response of nucleon 
matter to electromagnetic probes including correlated particle-hole rescattering have 
been relatively successful. In these works, as well as here, the three- and higher-body effective 
interactions are neglected. 

In the present work we use the static pair correlation operator: 


F., = E . 

p=l,6 


Of- — 1, Tj ■ Tj, (Tj • CT,-, Tj • Tj(Ti ■ (Tj, Sij, Ti ■ TjS., 






(9) 

( 10 ) 


5 















In place of the p = 1, 6 superscripts we often use the letters c, r, a, ar, t and tr to denote 
the radial functions associated with these operators. For example, 


fP=li6/ 'i _ fC fT ra rar rt rtr 

J \ ‘ij) — J ij'^ J ij’’ J ij’’ ^ ij ’ U ’ J ij * 


( 11 ) 


The Fij is obtained by minimizing the energy of symmetric nuclear matter at density p = 
Pn + Pp using hypernetted and operator chain summation methods [^, The results of 


the latest variational calculations are briefly summarized in section VI for completeness. 
The Argonne vl8 two-nucleon and Urbana IX three-nucleon [^] interactions are used in 
these nuclear matter calculations, in studies of weak interactions of light nuclei [IT, Q. and 
in the present work. However, improved models of Vijk are now available |^. The variational 
calculations of nucleon matter also include two spin-orbit terms in the Fij which are omitted 
here for simplicity. The variationally optimized Fij can depend upon the proton fraction Xp. 
However, this dependence seems to be relatively weak. The effective interaction obtained 
from the Fij in symmetric nuclear matter gives a fair description of the spin susceptibility 
of pure neutron matter. 

Matrix elements of operators between CB states are generally calculated using cluster 
expansions 0. We begin with the simplest, lowest order two-nucleon cluster approximation 
to study the general properties of the weak one-body effective operators and of the two-body 
interactions in CB for nucleon matter at densities p = 0.08, 0.16 and 0.24 fm“^ and for proton 
fraction Xp = Pp/p = 0.2, 0.3, 0.4 and 0.5. In this density range the contributions of clusters 
with > 3 nucleons to the energy of symmetric nuclear matter increases from 10 to 30% of that 
of the 2-body [^; thus the present results have only qualitative signihcance. We study the 
density, proton fraction and momentum dependence of the operators and the interactions. 

Due to correlations and weak pair currents, the effective weak current operators have 
2- and many-body terms in addition to the leading one-body term we consider here. The 
lowest order (in cluster expansion) effective one-body F, GT and neutral current operators 
are calculated and their results are presented in sections H to V. As expected the one-body 
CB matrix elements are smaller in magnitude than those in FG states. The dominant 
term responsible for the quenching arises from pion exchange interactions which change the 
isospins and spins of the nucleons. In the FG wave function, a nucleon in the single particle 
state for example, is a spin | neutron with unit probability. This probability is 

reduced in the CB state by the spin-isospin correlation operators acting on the FG state. 
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In contrast, the spin-isospin independent spatial correlations induced by the repulsive core 
in the two-nucleon interaction increase the magnitude of F, GT and NA matrix elements; 
however, they quench the neutron NV. The CB matrix elements of the charge current oper¬ 
ators are found to have a rather small dependence on the matter density and Xp within the 
range considered. They depend primarily on the momentum transfer q, and only slightly 
on the initial or hnal nucleon momentum. In addition to these, the neutral current matrix 
elements also depend upon Xp. The proton NV matrix element is an exception; it has large 
cancellations and depends on all of the relevant variables. 

The squares of the F and GT matrix elements in GB states are ~ 0.8 and 0.75 times those 
in FG states at small values of q. Thus the present 0th order (in GB Hj) 2-body cluster 
calculation predicts a quenching of low energy weak transitions in nuclei and nucleon matter 
by 20 to 25 %. It is likely that higher order effects will further reduce the matrix elements 
and increase the quenching. For example, the occupation probability of states with momenta 
;< /cf is ~ 0.9 in GB states, and it decreases to ~ 0.8 on including 2nd order Hj corrections 
3^ . In order to obtain quantitative results it will be necessary to include contributions 


of > 3-body clusters to the GB matrix elements neglected in this initial study. This has 
been done for symmetric nuclear matter with operator chain summation techniques; 
however, they are difficult to use in matter with Xp ^ 0.5. Three-body cluster contributions 
in asymmetric matter can now be calculated using the recently developed matrix methods 


The results for the GB two-nucleon interaction are presented in Sect. It gives a 

fair description of the spin-isospin susceptibilities of nucleon matter used to determine the 
effective interactions in the Landau-Migdal scheme |0] . It also has the typical features of the 


effective interactions used in existing calculations of weak interactions in nucleon matter [18 


If we assume that the calculations with effective interactions are implicitly using GB states, 
then their results should be reduced by a factor of ~ 0.75 to take into account the quenching 
of the F and GT matrix elements by short range correlations. Attempts to calculate the 
weak interaction rates in nucleon matter with the effective operators and GB interaction 
presented here are in progress. 
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II. CORRELATED BASIS FERMI MATRIX ELEMENT 


Let |/) and \F) denote the normalized correlated states obtained by operating on the FG 


states |<F/) and by the correlation operator SUFij. The CB Fermi matrix elements are 


given by: 





( 12 ) 


apart from the orthogonality corrections neglected here. The corresponding uncorre¬ 
lated, FG matrix element (FGME) is (<hir|Oi?|<F/). It is non-zero only when the occupation 
numbers of the states <F/ and <Fi? differ by only one nucleon, since is a one-body opera¬ 
tor. In contrast the GB matrix element (GBME) can be non-zero even when the occupation 
numbers of <F/ and differ by more than one nucleon. However, here we consider only 
the dominant “one-body” GBME in which they differ by only one nucleon. We dehne the 
quenching factor, rj, as the ratio of the square of these matrix element, |GBMEp/|FGMEp. 

We assume that |<h/) has full neutron and proton Fermi spheres with momenta kpn and 





(13) 


where kn < kpn and kp > kpp- In the absence of spin-orbit correlations, the Fermi matrix 
elements are non-zero only when the spin state Xn = Xp- The FGME=1 when kp — k„ = q. 
These conditions are also necessary for the GBME to be nonzero; however, its value can 
depend upon the matter density, proton fraction and the magnitudes kn, kp and q. 

The cluster expansion of the GBME is obtained by replacing the correlation operators 
Fij by 1 -|- {Fij — 1) 1^ and expanding the numerator and the denominator in powers of 
{Fij — 1). It is convenient to use the <Ff, containing only a product of single-particle wave 
functions in which nucleons i are in plane wave states with momentum k* and spin-isospin 
Xri'i), in place of the antisymmetric <F/ and use the antisymmetric This is equivalent to 
retaining the antisymmetric and <h p and has the advantage that we can associate nucleon 
numbers with the state labels in . The nucleon in the state k„y;„ of is labeled “a” for 
active; in uncorrelated states only a participates in the transition. All of the other nucleons 
in the Fermi spheres are denoted by j. 

The cluster expansion of the GBME is represented by diagrams as shown in Fig. |I[ The 
terms in the expansion are labeled with F.n.x.y, where F stands for Fermi, n is the order 
of the {Fij — 1) correlations, x = d,e for direct and exchange terms, and y = a, j denoting 
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the nucleon on which the weak interaction operates. The dots in these diagrams denote 
nucleons, a thin line specifying the states occupied by the nucleon in and passes 
through each dot. The nucleons a and j occupy states k„ and kj in the , therefore lines 
labeled k„ and k^ originate from them in all diagrams. Their termination depends upon the 
exchange pattern, since $^7 is antisymmetric. In direct terms the state line kj emerges and 
ends in the dot j because the state of nucleon j is unchanged. The line with the two labels 
k„ and kp denotes the weak transition. In direct diagrams it begins and ends in the dot 
a. In diagrams in which a and j are exchanged, the transition line begins at a and ends in 
j, while the state line kj begins from j and ends in a. The state and transition lines must 
form closed loops in all diagrams. The dashed line attached to nucleon i = a, or j shows 
the Fermi operator Opli) = The {Fij — 1) correlations are indicated by wavy lines. 

We sum over the spin-isospin states XrU) of the nucleon j, while those of a, Xn and Xp, are 
specihed by <hi7 (Eq. (0)). 

The equations for F.n.x.y are given below in the two-body cluster approximation in 
which n <2. They show that the F.n.x.y are independent of q, kn and kp when x,y = d, a; 
they depend only on q when x,y = d,j; and only on kn and kp in exchange diagrams 
{x = e). We also give a simple explanation of the important F.2.d.a term responsible for 
much of the quenching. The standard 2nd order perturbation theory calculation of the direct 
contributions to the Fermi matrix element is reviewed in Appendix A. One can easily identify 
the analogues of F.n.d.y in that familiar theory and obtain relations between the present 
approach and that of Arima and coworkers JT^. The perturbation theory assumes that the 
forces are weak, but in reality we cannot expand in powers of the strong, bare 2-nucleon 
interaction However, we hope that standard perturbation theory can be used in CB with 
the effective operators and interactions described here, as mentioned in the introduction. 

The leading 0th order term is given by: 


F.O.d.a = FGME 


d^r ^'’^■’■(Xp(a)|T+(a)|x„(a)) = 1 . 


(14) 


The momentum conserving delta function (5^(kp — k„ — q) and the Xn = Xp spin constraint 
are implied here as well as in all terms of the expansion given below. There are no other 
0th order terms. 
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The 1st order direct term with Opij 7^ a) is given by: 

F.l.d.j = ^ f d^raj {xp{a)xT{j)\{rj , {Faj - l)}\Xn{a)xT{j)) 

j 

= p j d^r 2 /"(r) . (15) 

All spin dependent terms in Faj give zero contribution on summing over the spin states of 
nucleons j, and the factor of 2 in the above equation comes from: 


{■rj , Tj ■ T„} = 2r+ . 


( 16 ) 


From now on the aj subscripts on r and F will be dropped for brevity, and the r dependence 
of the /^’s will be implicit. 

The contribution of F.l.e.j is given by: 


F.l.e.j = 




3r ‘'"^■’'{Xp(a)Xr(j)|eai{T+ , (F - l)}|x„(a)x^(j)) 


= - dh [pn in{r){r - 1 + 3 /") + Pp ip{r){r + Sf 


(17) 


where Cij is the spin-isospin exchange operator: 


Cij = - j(l + r, ■ Tj)(l + <Ti ■ <T,) , 


( 18 ) 


and the Slater functions (A^ = n, p) are: 
2 f d^k 


iN{r) = 


Pn j (27r)3 


9{kFN — k) = ?)[sin{kFNf’) — kFN^ cos{kFNr)]/{kFNfy ■ (19) 


The algebra of the operators OF described in Ref. is very useful in evaluating these 


contributions. 

The 2-body terms with Of{(i) have contributions from the numerator of the matrix 


element, Eq. (^), as well as normalization corrections introduced through the expansion 
of the denominator. We denote these by F.l.x.a.N and F.l.x.a.D respectively. In Fig. |l| 
the denominator contributions are shown as products of two diagrams. The 1st order direct 
terms with Of{(i) cancel: 


F.l.d.a = F.l.d.a.N + F.l.d.a.D = 0 , 


( 20 ) 
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while for the exchange terms we obtain: 


F.l.e.a = F.l.e.a.N + F.l.e.a.D (21) 

F.l.e.a.N = ^ f d^r , (F - l)}|xn(a)Xr(j)) 

j 

= -Jd^r F^^-^[ppip{r - 1 + 3/") + Pninir + ^D] , (22) 


F.l.e.a.D = / d^r(e ‘'"^'’'(Xn(a)Xr(j)|eai(D - l)|xn(a)Xr(j)) 

j 

+ {xp{a)xT{j)\eaj{F - l)|xp(a)Xr(j))) 

= d^r{[r-i + sr + r + ^r] l + ip] 

+2[r + + F*^- Vn4]) . (23) 

For calcnlating the 2nd order terms, it is convenient to dehne: 


F = i + f° + FV,-t,- , 

(24) 

F° = /= - 1 + rCTa ■ CT, + fSaj , 

(25) 

= r + + f-Sp ,,. 

(26) 

Only the spin independent parts of the prodncts of the above F° 

and F^ contribnte to the 

second order diagrams. These are called the C— parts in Ref. 31 

. We dehne: 

cy = C[F‘F^] , 

(27) 

= ^[(l + a, ■ (t^)F^F^] . 

(28) 


The expressions for and in terms of the correlation fnnctions, /^, are given in 
Appendix p. 

There is no contribntion from the denominator to the terms F.2.x.j. These are given by: 

F.2.d.j = ^ f d^r {xp{a)Xr{j)\{F - l)Tj{F - l)\xn{a)Xr{j)) 
j 

= p J2[cy + Cf] , (29) 


F.2.e.j 


j dV ‘"^^■'■{xp(a)x^(j)|eai(F - 1 )t+(F - l)|x„(a)x^(j)) 
j 

_i J rfVe“‘..-[p„(.„(C«»-Ci') + 2p„(,(C“+Cf)] . (30) 
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The sum: 


F.2.d.a = F.2.d.a.N + F.2.d.a.D = 

E J d^r{xp{a)xr{j)\{F-1 )t+{F-1 )(F - l)2}|x„(a)xr(j)) 

= p [ d^-AC}^) . 


(31) 


Note that only the F^Ta ■ tj, which does not commute with the operator, contributes to 
this sum. 

The results presented in the next subsection show that the above term gives the largest 
contribution to the quenching of the Fermi matrix element in matter. This term simply 
takes into account the probability for nucleon a to be a neutron in the initial and a proton 
in the final state. In the uncorrelated product state, |*hf), nucleon a is n but in the 
correlated product state, 5nFy|<Ff), it can be in other nucleon states. We refer to nucleon 
a in the correlated state as a “quasi-nucleon”. The probability that it is a neutron is given 
by: 

(«>,|[5nF„] i(l-r') |5nF.,]|4.f) 


P/(a = n) = 


(32) 


{<I.,|[5nf’y]2|<l.f) 

We use the cluster expansion to calculate this probability. The 0th order, one-body term is 
unit, and the two-body 2nd order direct term is: 

- IE /<i’>'{x'„(o)\v(j)|(F-l)r.'(f-l)-lK. {F-lf}\Xn{a)XrU)} 


^ Pp d r(-4Cj ) . 


(33) 


The two-body 1st order direct terms cancel as in Eq. (pOl). Neglecting the exchange terms 
we obtain the direct part: 


Pi{a = n, d) = 1 + pp / Pr(—4(7]^) . 


(34) 


In a similar way, the direct part of the probability for the active quasi-nucleon, a, to be a 
proton in the hnal state is given by: 


Pp[a ^ p, d) ^ 1+ p„ / d'*r(-4C] ) 


(35) 


Hence 


1 -I- F.2.d.a = Pi{a = n, d) Ppicb = p, d) , 


(36) 
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neglecting the terms of order [C^Y- 

The probabilities for the active quasi-nucleon to be in the initial spin isospin states t, J. 
n,p have been calculated keeping only the direct terms, at the three densities for Xp = 0.5. 
These are given in Table |. In one-body Fermi transitions these are also the probabilities 
for the active quasi-nucleon to be a spin I p,n in the final state. 

The 2nd order exchange term: 

F.2.e.a = F.2.e.a.N + F.2.e.a.D , (37) 

has contributions from both F^ and They are given by: 

F.2.e.a.N = ^ f d^r F^'"p~'^^'>'YXp{a)XT{j)\eaj{F - 1 )t+{F - l)\xn{a)xT{j)) 
j 

= -\j [pp ipiCY - CY) + 2p. in{CY + , (38) 

F.2.e.a.D = / d^r(e-"^^^-'""'>'Yxn{a)xT{j)\eaj{F - lY\Xn{a)Xr{j)) 

+ e-^^^^-^^^-YXp{a)Xr{j)\eaYF - lY\Xp{a)Xr{j))) 

= \J d^r(^{-ACY+AClY{F^--^Ppip + F^^'’^Pnin) 

+{CY + CY + 2CY){F^--^p^in + F^-'^Ppip)') . (39) 


A. Results for Fermi Matrix Element 

The Fermi matrix elements have been calculated using correlation functions obtained in 
Ref. by minimizing the energy of symmetric nuclear matter using the Argonne-vl8 and 
Urbana IX 2- and 3-nucleon interactions. In Fig. ^ we present the results for pp, the square 
of the Fermi CBME (Eq. (|T^), for kn = kpn and kp = kpp. 

When Xp < 0.5 the total isospin Tj of the state |J) is [N — 2’)/2, while that of |F) is 
{N — Z)/2 — 1. In the case of symmetric nuclear matter the T/ = 0, while Tp = 1. Thus the 
calculated matrix elements are between states with AT = 1. The Fermi matrix elements 
for g = 0, between isobaric analogue states having the same T and Tzf — Tzi ± 1 are given 
by (T =F Tzi)(T ± Tzi + 1) in both correlated and uncorrelated states. We will not discuss 
AT = 0 Fermi ME in this paper. 
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The variation of tjf with proton fraction is less than 3% at all densities calculated. How¬ 
ever, the proton fraction limits the allowed values of q through the momentum conservation 
relation: q = kp — k„. The variation with total density is also small within the considered 
range. This suggests that we can approximate |CBMEp by a function of p and q. In the 
small q region, q < 0.5 fm“^, it can be well represented by the quadratic: 

Vf = r]F{q ^ 0) + apip) q^ ■ (40) 


We have £t the calculated values for symmetric nuclear matter and the results are given in 
Table |TI|. 

Fig. ^ shows the contributions of each term in the cluster expansion of the Fermi matrix 
element in matter at density po and Xp = 0.5. The F.n.x.a and F.n.e.j terms give contribu¬ 
tions that are independent of q as can be seen from Eqs. (0), (HD, (H), ®, (H), m 


and (^). The dominant contribution to the quenching of the Fermi CBME comes from 
F.2.d.a] |1 -|- F.2.d.a\^ = 0.7 is shown by the dotted line in Fig. ^ As discussed in the 
previous subsection this result can be interpreted in terms of the probabilities for the active 
quasi-nucleon a to be a neutron in the initial and a proton in the final CB states. 

The exchange terms, F.n.e.a, contribute an additional ~ 0.1 to the g-independent 
quenching; | '^n^F.n.x.al'^ = 0.61 is shown by the double dash-dot line. This additional 
quenching is mostly canceled by the F.n.e.j terms, as shown by the dash-double dot line; 

I + F.n.e.j\‘^ = 0.71. 

The F.n.d.j terms, given by Eqs. (|T^ and (HD) introduce the g-dependence. Of these, 
the 2nd order F.2.d.j is dominant as can be seen from the dashed line, which includes only 
F.l.d.j and all of the g-independent terms. The full line gives the square of the total matrix 
element including F.2.d.j. 

The contributions of the various correlations to the CBME are shown in Fig. ^ The 


1st and 2nd order terms are dominated by the /'^’^(r* 


(TjTi 


and 


correlations induced mainly be the OPEP. After setting = 0 the | '^nxF-n.x.a\‘^ 

becomes essentially 1 as shown by the dotted line in Fig. The full CBME exceeds unity 
in this case (see the dashed line) via the contributions of Z^ — 1 correlations to F.n.x.j. The 
dash-dot line shows rjp obtained by further setting Z'^ = 1- It is fairly close to one showing 
that the and Z* correlations have small effects. 

The Fermi CBME, calculated in the two-body cluster approximation does not depend 
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significantly on the magnitndes of initial and final nncleon momenta. The dependence on 
kpn — kn and kp — kpp is illustrated in Fig. It shows rjp for p = po, Xp = 0.5 and 0.3, 
kn = 1, 0.75, 0.5 kpn and kp = 1, 1.25, 1.5 kpp as a function of q. The results for the 18 
possible combinations of Xp, k^ and kp values differ by less than 0.03. 

III. CORRELATED BASIS GAMOW-TELLER MATRIX ELEMENT 

The procedure for the calculation of the GT matrix element is similar to that discussed 
previously in Section We therefore discuss only the differences and give the hnal ex¬ 
pressions. The operator, Ogt, is an axial vector and it is convenient to express its matrix 
element using the following two axial vectors: 

(^) = (Xp(a)l o'(a) 'r+(a) \Xn{,a)) (41) 

and 

(At) = 3 Taj i^a) ■ raj “ {^a) , (42) 

obtained from the tensor correlations between nucleons a and j. Note that (A^) depends 

upon Taj. 

We assume that Xn in Eq- (0) is spin up and sum the square of the GT matrix element 
for the two hnal states with Xp =1,1 denoted by \F |) and \F 1). In FG we get contributions 
only via the operator ax = o'aT^] only a^ia) contributes to the FGME with Xp =T) while 
ax{a) and o'y(a) give the GT FGME for Xp =i- However, in GB the A^ induces transitions 
that are forbidden in FG states. 

The terms in the cluster expansion of the GT GBME are denoted by GT.n.x.y as in 
the last section. The ratio qa of the axial to vector coupling constants is omitted from the 
GT.n.x.y for brevity. We obtain: 


GT.O.d.a = (cTq) 


(43) 


GT.l.d.j 

GT.l.e.j 


p jd^re-^^-^ 

+Pii ip — 


2 (r^ ( 5 v) + /‘" (A,)) 

{firXr+ipn i^a) 
i + r + 2/") (Sv) + (/' 



(44) 

(45) 
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GT.l.e.a.N = - ^ 4(r + 3/"^ (^) 

+p, 4 [ir -i + r + 2 r) (^) + if - n ( a ,) 

GT.l.e.a.D = {d^)F.l.e.a.D 


(46) 

(47) 


GT.2.d.j 

GT.2.e.j 


2p dhe-“>GF^,y + F^p) (Sv) + (Fiy + F”^) (A,) 


(48) 


-- d\e‘^-'' 


pr 




00,A 


+«5’^+^n (A.) 


pll,A _ plO,A 
01,A\ 


(F”^" - F^’ - F^f + FiY) (S;) 

+ GG (A,)! + 2ft 4 [k;’" + Gf) ('^) 


(49) 


GT.2.d.a 

GT.2.e.a.N 


GT.2.e.a.D 


P d^r{Fl 


00,(7 

a 


Fpf - cr - 3cy) (a,) 


(50) 


— Fre 


ikri-r 


/ r00,O- _ pll,(T I rlO,(T _ p01,(T\ /~\ 
\‘^e,a e,a ' e,a e,a ) 


+ ( 4 “/ - Go!' + G/' - (A.) + 2p„ 4 (FlY + F“'”) (<4) 


+ (P,a'" + P"/) (Ai) 

{aY)F.2.e.a.D 


(51) 

(52) 


The coefficients (p = a,j) and F^p are defined as the cr^ and At parts of the operator 
F^a-yF-^: 

F^a-yF-^ = F^y'^CTa + F^pAt + terms linear in cr^ . (53) 

The remaining parts linear in cr^ do not contribute after summing over (jj. The F^y'°^ and 
Fl^’^ are the corresponding parts of the operator (1 + cr^ ■ (Tj)F^cryF'^, and the expressions 
for 44'" and 44 ^ are given in Appendix p. 

As in the Fermi case, the 2nd order direct diagrams, GT.2.d.a can be interpreted in terms 
of quasi-nucleon probabilities. The GT.2.d.a has contributions from the (4) only. When 
the final proton has spin | only the 4 term contributes. We consider this simple case for 
illustration. In this case, (1 -|- GT.2.d.a) represents the probability that the active quasi¬ 
nucleon has 44 ~ the initial state and -|-1 in the final state. In FG states these are 

unit probabilities. We use the cluster expansion to calculate them in CB states. The 0th 
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order terms = 1, and the two-body 2nd order direct terms are given by: 


= I (F^T - - CT - ^CY) 

1 


T (Pp -Pn)^ I dh + 2F];/ - 2CY + 2CY) , 


(54) 


where the upper and lower signs correspond to the initial and hnal states respectively. The 
1st order direct terms cancel as in Eq. (|20|) . Neglecting the exchange terms and those of 
order and we obtain; 


Pi{« = 


-1, d) P^a: r: = 1, d) = 1 + p y dV - FJ;/ - Cf - 3CY) 

= 1 + [GT.2.d.a]^ . (55) 


The cr^ rf = —1 probability is the sum of the n ] and p | probabilities listed in Table 


Results of Gamow-Teller Matrix Element 


The tensor correlations lead to a dependence of the GT CBME on the direction of the 
spin quantization axis through the (A^) terms. We therefore do not discuss the CBME for 
spin up and down hnal states individually. The sum of |CBMEp over the hnal two spin 
states determines the transition rates and is independent of the chosen axis. This sum equals 
3 for FGME. In the following we report results for: 


>)G7- = 3(|{rT |OGT|/>P + l(ri |Ogt|/>R 


( 66 ) 


The rjcT has been calculated using the correlation functions as described in Section |11 A 
and the results for = kp^ are plotted in Fig. ^ As in the Fermi case, the variation of 
Pgt due to changes in proton fraction is less than 3%, but it has more q dependence. The 
quadratic ht (Eq. (^0|)) is still valid up to g ~ 0.5 fm“^, and its parameters are given in 
Table |TI|. 

Fig. 0 illustrates the relative contributions of the various terms to Pgt- As in the Fermi 
case, the main quenching comes from the GT.2.d.a term; approximating the |CBMEp by 
|1 -|- GT.2.d.a\‘^ gives pgt = 0.79 (dotted line). It decreases to 0.72 on adding the GT.n.e.a 
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terms (double dash-dot line). The double dot-dash line shows the result after including 
GT.n.e.j terms which reduce the quenching. 

The main g-dependence comes from the 1st order GT.l.d.j term; results obtained after 
adding this term are shown by the dashed line. The GT.2.d.j term also contributes to the 
g-dependence (full line gives the total tjgt)- 

The dash-dot and the double dot-dash lines have a barely visible g-dependence coming 
from the GT.n.e.y terms. In the Fermi case these exchange terms depend only on kn and 
kp] however, in the GT case they introduce a dependence of rjcT on the angle between k„ 
and kp. This appears as a g-dependence, but it is very small (< 0.002). 

The relative contributions of various correlations to rjcT cU'e shown in Fig. |^. The dashed 
and the dash-dot lines show results obtained after setting = 0, and in addition 

f^ = l respectively. The central correlations contribute mostly via the GT.n.e.j terms; the 
dotted line close to rjcT = 1 is obtained by setting = /^"^ = 0 and including only the 
GT.n.x.a terms. 

The dependence of r]GT on kn and kp is shown in Fig. It is small < 0.03 as for r]F- 

IV. CORRELATED BASIS NEUTRAL-VECTOR MATRIX ELEMENT 

In “one-body” NV transitions the hnal state is: 

|$i?) — |*h/) , (57) 

where fc* < kp^i and kf > kpM- The NV matrix element is nonzero only when the initial 
and hnal spin-isospin states, Xn Xv same. 

The terms in the cluster expansion of the NV CBME are denoted by 

NV.n.x.y = —sii)?9w NV.n.x.y.l + -(1 — 2siT?6w) NV.n.x.y.z , (58) 

where n, x and y are dehned in section y and 1 and z respectively represent contributions 
of and For the NV.n.x.y.l terms we obtain: 

NV.D.d.a.l = 1 (59) 
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(60) 


NV.l.d.j.l = /'<iVe-<‘-'2((/“-l)p + r(p,-p„){r'» 

W.l.e.j.l = -i j<ere'^"{{r-l + 3r + 3r + 9r){p^e, + PnQ 

+(/' -1 + sr - r - 3r)(pp(p - pntn)[K)] (ei) 


NV.l.d.a.l = 0 (62) 

NV.l.e.a.N.l = -I jdh [{f - I + ?,r + + p^L) 

+(r -1 + 3r - r - srniPp^ - PnC)(r,^)} (63) 

NV.l.e.a.D.l = i Jd^r + e^'^n[{r - 1 + Sr + Sf + 9r){Pp^P + PnL) 

+(r -1 + 3r - r - 3r){pp^ - PnQir:)] (64) 

NV.2.d.j.l = |dVe-q-[(Cr + 3C]i)p + 2(C]°-C]i)(p,-p„)(r:)] (65) 

NV.2.e.j.l = y'd3re**^--'{(C'“ + 6C'f-3 C']')(Pp4 + P„C) 

+ (C“ - 2Cl° + 5Cl^){p,i, - PnDiTl)] (66) 


NV.2.d.a.l = 0 (67) 

NV.2.e.a.N.l = ~ Jd^r + QCl^ - 3Cl^){ppip + p^L) 

+ (C“ - 2Cf + 5Cl^){ppip - Pnin){T:)} (68) 

NV.2.e.a.D.l = I [(^oo ^ g^oi _ 3Cl^)^ppip + pnL) 

+ (C“ - 2C°i + 5Cl^){ppip - PnL){r:)] (69) 

where (t^) = {xUa)\ra\xM)■ 

The NV.n.d.a.l terms are zero for n > 0 because commutes with the static correla¬ 
tion operators. Also note that the exchange, NV.n.e.a.l, terms are zero when |kj| = |k/|. 
The NV.n.x.y.z terms are given by: 


NV.O.d.a.z = {t^) (70) 

NV.l.d.j.z = jdhe-^^-^ 2 ((r - l)(Pp - Pn) + /V(r:)) (71) 

Nv.i.e.j.z = y'rfVe**^-^|(r-i + 3r + r + 3r")(Pp4-p„c) 

+ (/'^ ~ 1 + 3/°^ + /^ + 3/‘^'^)(pp£p -I- PnC)(Ta)| (72) 


19 



(73) 


NV.l.d.a.z = 0 

NV.l.e.a.N.z = jd\ [{r - I + ?,r + f + ?,r){pptp - 

+(r -1 + 3r + r + ^DiPp^+Pn c)(r:)} (74) 

NV.l.e.a.D.z = {T^^)NV.l.e.a.D.l (75) 


NV.2.d.j.z 

NV.2.e.j.z 


1 


(0S° - c“)(p, - p„) + 2(C‘» + C“)p{r^) 

(C“ + 2C“ + C“)(p,|-,-p„4) 

+ (C“ + 2C“ + Ci‘)(p, + p„ ('„)(t;) 


■j jrfre*^" 


(76) 


(77) 


NV.2.d.a.z = y‘rfV(4C]i(p,-p„)-40(r:)) (78) 

NV.2.e.a.N.z = y’dVe**^^"-|(C'“ + 2^]° + C']^)(pp4 - p„4) 

+ ((7°° + 2C]° + C']^)(pp + Pn ^n){Ta) | (79) 

NV.2.e.a.D.z = {T^)NV.2.e.a.D.l (80) 

In symmetric nuclear matter the matrix elements of are related to those of r^. In this 
case the NV.n.x.y.z = F.n.x.y. However, when Xp < 0.5 the NV matrix elements have 
additional terms dependent upon pn — Pp, or equivalently Xp. 


A. Results of Neutral Vector Matrix Element 

In uncorrelated FG states, the neutral vector matrix element is: 

-sin^Ow+ ^{l-2sm'^9w){Ta) = -0.2314 ± 0.2686 (81) 

for proton and neutron particle-hole pairs respectively. The above two terms nearly cancel 
for uncorrelated protons. The correlations influence each operator differently and the hnal 
CB result depends sensitively on ki,kf,p and Xp. The strong dependence of the proton 
NV matrix element on p and Xp is shown in Fig. where we have plotted the proton 
particle-hole NV CBME scaled by 0.0372, the FGME. Note that the value of the GBME 
(not \CBME\^) is shown in this hgure. At low densities the hrst term dominates, and the 
GBME is negative; however, at higher densities the second term becomes larger, and the 
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matrix element becomes positive. At p ~ po the cancellation of the two terms is almost 
exact, and the proton NV CBME is very small. Fortunately, in this case the FGME is small 
and the CBME is of the same order in the considered density range. Thus, the coupling 
of the proton NV current is not likely to have a signihcant contribution to the i/-nucleus 
interaction. 

Figure 0 shows the density and Xp dependence of ptvv for neutron particle-hole pair 
excitations. At p = ^po the correlations increase the contribution of the hrst term and 
decrease that of the second term in Eq. (ED by a similar magnitude. Therefore at small 
p and q the NV neutron CBME ~ FCME. However, at higher densities it is quenched. 
As mentioned earlier these matrix elements have a signihcant Xp dependence absent in the 
charge current matrix elements. 


Fig. |T^ shows the contributions of the various correlations to the NV neutron CBME. 
The CBME is inhuenced by contributions of the — 1 correlations to NV.n.x.j.l and those 
of the ■ cTjTi ■ Tj and f^^{r-ij)SijTi ■ Tj correlations to NV.n.x.y.z terms. The 

results obtained after setting = 0, and in addition /'^ = 1 are shown by dashed 


and dash-dot lines in Fig. 12. 

The Neutral Vector CBME for a neutron particle-hole pair does not depend signihcantly 
on the magnitudes of the initial and hnal nucleon momenta. Variation of ki from 0.5 to 1 
and of kf from 1 to 1.5 kpn changes pAry by less than 3 %. 


V. CORRELATED BASIS NEUTRAL-AXIAL-VECTOR MATRIX ELEMENT 


The operator. On a is an axial vector and it is convenient to express its matrix element 
using the following two axial vectors, similar to those used for the Camow-Teller CBME 
(Sect. Ill); 

(<^a) = (Xw(a)l cr(a)lx)v(«)) (82) 


and 


(At) = 3 Yaj {(Ta) ' Taj - {(Ta) 


(83) 


We assume that x^n i’^ ^0- (1^) is spin up and calculate the sum of the square of the NA 
matrix element for the two final states with Xn i both N = n and p. The terms in 
the cluster expansion of the NA CBME are denoted by NA.n.x.y as in Section |1|, and the 
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factor qa is omitted for brevity. We obtain: 


iVAO.d.a = (84) 

NA.l.d.j = ^ 2[{pp-pn){r{(Ta) + 

+p{r:){r {a,) + r {A,))] (SS) 

A^Al.e.j = jd\e^^-^[{{r -1 + r + r -?,r){aa) + {f + ?,n{A^)){pp^p-pnL) 

+((r -1 + r + r + ^r){^a) + (/* - n{At)){ppip+ p„c)(a^)} (se) 


NA.l.d.a 
N AA.e.a.N 


= 0 




(87) 

i + r + r- 3r){<T.) + (/' + 3/‘-)(A,))(pp«, - p„« 


+((r -1 + r + r + 5 /"){<t.) + (/• - /■-)(A,))(p,«p + p„<„)(r„")} (ss) 

TVyf.l.e.a.D = -(cra)(r^)iVKl.e.a.Zi).l (89) 


NAA.d.j = - (a,) + (F°f - (A,)J (p, - p„) 

+ (A" + db + 2Fiy) (o-J + (F]"’-" + + 2F“--<) (A,)] p(r') (90) 

1 


tttOO.Ct I o ITilO.fJ 


^01,0- 

ej 


NA.2.e.j = 

+Kr+3i^r-^i^+^ii(At 

- Fi°’" + 3F^y + F^y) {aa) 


V — n P \ 

p ''p Hn ) 


ej 


ej 


+(d; 


ej 


+ 3F“--" + p}f) (A,)) (p, 4 + p„ 4) (4) 


(91) 


NA.2.d.a = 


'd^r 


(FT + fff + 2db - 2(CS‘ - C]')) (^„)(p, - p„) 


NA.2.e.a.N = — /rfVe“' 


+(d,r - dr - c“ - 3cr) (-^.(pw) 

1 


^kf-r 


p00,C7 plO,C7 Q pOljCT 

■^e,i “ -^ej + 


ilO,a 


HK 


00,A 


F!^P^ + 3Fy^ + Fy^){At 


e,J 


F^^ + 3Fy 


10,(7 


ej 
701,(7 

e.i 


+di )K> 

Pn 


P 


+dj )K) 

+(F“y+spp - F)',y+F.‘y) ( a,) } (p^ 4 +p„ o ( 4) 


N A.2.e.a.D = -{aa){T^)NV.2.e.a.DA 


(92) 


(93) 

(94) 
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A. Results of Neutral Axial-Vector Matrix Element 


We discuss only the sum of the |CBMEp over the two hnal spin states because it is 
independent of the chosen spin quantization axis. This sum equals 3/4 for FGME. In the 
following we provide results for: 


VNA = MF T |Om|/)P + \{Fl |0„4|/)n 


(95) 


The rjjsfA for neutron and proton particle-hole pairs are plotted in Figures O and P 


re¬ 


spectively for the considered density and proton fraction values. In these matrix elements 
ki = kf = kpN- 

The charge-changing and neutral axial vector operators (Ogt and On a), appropriately 
scaled, can be interpreted as the three components of an isospin vector operator. In sym¬ 
metric nuclear matter the expectation values of these three components are equal as one can 


not quantify the isospin axis. The stars in Figures and |T^ are results obtained for rjcr 
for symmetric nuclear matter with equivalent initial and hnal momenta and densities. They 
are identical to those obtained for tjna for both proton and neutron particle-hole pairs. 

Unlike the results for the GT GBME, there is a noticeable dependence of t]na on the 
proton fraction at all densities considered. This {pp — Pn) dependence originates from the r/ 
in NA.n.x.j and NA.n.x.a terms. We can approximate the NA results obtained for Xp < 0.5 
by adding a density dependent term proportional to {pp — pn) to pna for symmetric nuclear 
matter. For small q this approximation is: 


Vna{p,Xp < 0.5) = Pna{p,Xp = 0.5) - Cn{p){Pp - Pn) (96) 

= Vgt{(1 = 0) + acTff — Cn{p){Pp — pn) , (97) 

where we have used tjna = Vgt at Xp = 0.5 and Eq. (^0|) . The values obtained for Cn{p) at 
the three densities considered are given in Table |T|. 

The correlation dependence and initial and hnal momenta dependence studied for pcT 
are applicable here and will not be discussed further. 


VI. CORRELATED BASIS INTERACTION 

The expectation value of 77 — Tfg (X), where Tpci^) is the kinetic energy of the Fermi- 
gas state <hx, is expanded to calculate the energy of the correlated state |X). It is given 
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by: 


(X\H\X) = + T..(x), (98) 


(^>x|[5nF,,]2|$x) 

J±_ 

2m 


TpciX) = Y1 ^ 

all i occupied in 

Since <l>x is an eigenstate of the kinetic energy operator T = Yli ~X‘f/2m, with eigenvalue 
Tpg{X), it is not necessary to expand the FG kinetic energy. The {H — Tpg{X))\X) does 
not contain terms with V? operating on |$;!(:)• Including only two-body clusters we obtain: 

i(V2fy) - -I 

m m 

where \ij) = The gradient operates on the relative coordinate, and 

the sum i < j is over states occupied in <Fx. The effective correlated basis two-nucleon 
interaction (CBI) is given by (see Eq. (8)): 


{X\H\X) = Tpg{X) + 


VijFij 


:vF,,) ■ V 


\^J) , (100) 


i<j 


..CB 


= F 

1-1 


VijFij 


-(V%) - -(VFy; 

m m 


■ V 


( 101 ) 


-IJ * o 

in the 2-body cluster approximation. The energies of correlated states |X) are obtained by 
using this in 1st order with FG wave functions, as in the Hartree-Fock approxima¬ 
tion. 

The has a momentum dependence via the (VF^j) ■ V term which gives contributions 
to the matter energy via exchange terms in Eq. (|100|) . This contribution is much smaller 
than that of the momentum independent, static terms in vfF defined as: 


CBS 

o 


— Fij{yij V )Fij 


m 


( 102 ) 


In the present work we have considered only the static part of Fij as mentioned in the 
introduction. We therefore keep only the dominant, static part of the full Argonne The 
full Vij is first approximated by a Ug interaction chosen such that it equals the isoscalar part 
of the full interaction in all S and P waves as well as in the ^Di wave and its coupling 
to the ^Si- The difference between the full and the Ug interactions is small and treated 
perturbatively in the quantum Monte Garlo calculations [^. The v'g has terms with the six 
static operators, and two spin-orbit terms. The later two are omitted to obtain the 

static part of Argonne In this approximation the is a static operator having six 

terms with 


= 


p=l,6 


,CBS, 






(103) 
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The Landau-Migdal effective interactions used in studies of weak interactions in nuclei 
and nucleon matter are obtained from the spin-isospin susceptibilities of nucleon 
matter. We have therefore studied these susceptibilities with the and . The energy 


of nucleon matter with densities and can be expressed as: 

E{p, X, y, z) = Eo{p) -h Er{p)x^ + E^{p)y‘^ E„r{p)z^ , 
X = (Pnt + Pnl ~ Pp1 ~ Ppl)/ P 1 

y — (Pnt ~ Pni + Ppt ~ Ppi ) I P 1 

Z = (Pnt ~ Pni ~ Ppi + Ppi ) /P ■ 


(104) 

(105) 

(106) 
(107) 


The r, cr and err susceptibilities are proportional to E~l^^, and Eq{p) is the energy of 
symmetric nuclear matter with x = y = z = Note that the E^{pq) is the familiar 
symmetry energy in the liquid drop mass formula. In principle, the above expansion is valid 
at small values of x, y, and z] however, within the accuracy of available calculations it seems 
to be valid up to a; = 1 |^, ^ . 


We have calculated the Et-^„^„t{.p) using the obtained from the Eij at p = |, 1 and 
I po- The results obtained with the n*"'® are given by full lines in Fig. while those with 
the simpler by dashed lines. The momentum dependent part of gives rather small 
contributions which may be neglected in the first approximation. The has a density 
dependence due to that of E^j. However, it has very little effect on E„ and E^t-] the results 
obtained from the | , 1 , | po essentially overlap. The density dependence of has 

a small but noticeable effect on the symmetry energy Er{p). 


The stars on Fig. show the values of Et-{p) extracted from recent variational calcu¬ 
lations of symmetric nuclear matter (SNM) and pure neutron matter (PNM) with the 
Argonne vl8 and Urbana IX interactions, assuming that Eq. (|104|) is valid up to a; = 1 for 
y = z = 0. The two-body seems to provide a fair approximation to the E^. 

We also consider the spin susceptibility of PNM given by the inverse of E^^^{p) dehned 
as: 


E 


PNM, 


(p,p) = Eo™(p) + E™(p)p^ 


(108) 


The results obtained with the and are shown in Fig. along with those obtained 
from quantum Monte Carlo calculations [^] with the static part of Argonne vl8 and Urbana- 
IX interactions. The two-body using F),- of SNM gives fairly accurate values of . 
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Fig. 0 shows Eq{p) and Eq^^{p) calculated from the at the three values of p. 
The stars in this hgure give results of the recent variational calculations with the full 

Argonne vl8 and Urbana IX interactions. At low densities the two-body is not a bad 

approximation; however, the Eq{p) obtained from it does not show a minimum at po- The 
3-body interaction and cluster contributions are repulsive and are essential to obtain the 
minimum. 

The 2-body is more accurate in predicting the susceptibilities than the equation of 
state, Eq{p). This is partly because the contributions of Tpc and to the Er^aprip) add. 
The contribution of Tpc to the is shown in Fig. |^, it is about half of the total. 

For this reason even relatively simple estimates [^| of are not too different from the 

current state of the art |^. In contrast, in SNM the large negative cancels the Tpc to 

produce a relatively small binding energy. Therefore the many-body clusters are relatively 
more important in the calculation of Eq{p). 

The results of the recent SNM calculations, which provided the Eij used here, are sum¬ 
marized in Table |I^. The 1- and 2-body cluster contributions are calculated exactly. The 
calculation of the 3-body cluster contributions from the static part of Eij are also exact. 
However, the 3-body contributions from spin-orbit correlations and forces, the n>4-body 
contributions and the difference between the variational and the ground state energies are 
estimated. The empirical Eq{p) assumes po = 0-16 fm“^, Eq{pq) = —16 MeV and an incom¬ 
pressibility of 240 MeV. The difference between the calculated and the empirical values is 


likely to reduce when the more realistic Illinois Vijk is used in place of the Urbana-IX. 
However, a part of this difference is due to the approximations in the calculation. 

Next we consider the non-diagonal CBI. Let a Fermi-gas state |<Fi?) differ from I*!)/) in 
the occupation numbers of two single particle states: 

\^p) = ala\^ajai\^i) . 


The matrix element of H between the CB states is given by: 


{E\H\I) = 


(109) 


( 110 ) 


The numerator of this matrix element contains terms in which the kinetic energy operator 
acts on the $/. They give: 

(4>^|[5nF,,][5nF,,]T|$,) 




= TpGiI){E\I) 


( 111 ) 
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When the correlated states are orthogonalized this term is zero. Neglecting it the two-body 
cluster approximation of the above matrix element is obtained as: 


{F\H\I) = {mn\ 


.CBS 


m 


- (V'-(FV')F + F(VF)-V) \ij) 


= {mn\v^^\ij) . 


( 112 ) 


The V' operate to the left while V to the right. When the momentum dependent term is 
negligible, this matrix is just the Fourier transform of . Using the algebra of operators 
0^2^'^, and Eqs. 


and (|103|) we obtain: 


-r"= 


g,r,s,t=l,6 

Here we have used: 




g,s=l,6 


V^f(rtj)S,j = S.j + V^f‘(r,j)J , 


( 114 ) 


and the matrices are given in ref |^. The Fourier transforms of the are given in 


Figures to Note that Sij = Scr^q cr^q — cTi ■ ctj in momentum space. 

The effective is weaker than the bare v, particularly at large values of q, as shown in 
Figs. P!B| - PD| . Perturbative corrections typically involve a loop integration over the momentum 
transfer q with a phase-space factor. Hence in these hgures we compare q^Vp^^{q) with 

q%p{q). 


VII. CONCLUSIONS 

We have calculated the effect of short range correlations on nuclear weak interaction 
matrix elements. At low energies and small values of q the charge current, weak transition 
rates are quenched by ~ 20 to 25 % in the simplest 2-body cluster approximation in 0th 
order CB theory. This quenching is relatively independent of the density and proton fraction 
of nucleon matter as well as the momenta of nucleons in the | to | po range. However, it 
depends upon the momentum transfer q. 

The dominant part of the quenching is due to spin-isospin correlations induced by the 
OPEP in the bare interaction. The OPEP changes the isospin of nucleons. For example, 
in the n —>■ p weak transition between uncorrelated states the active nucleon is initially a 
neutron and hnally a proton with unit probability. In correlated states these probabilities 
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are less than unit, and they reduce the weak interaction matrix elements. In particular, 
for the Fermi case, most of the q independent reduction is given by the product of the 
probabilities for the active quasi-nucleon to be initially a neutron and hnally a proton. A 
similar interpretation is also applicable for the GT matrix elements. 

In contrast to charge current, neutral current matix elements have a signihcant depen¬ 
dence on the proton fraction. The neutron NV matrix element also depends on the total 
density, while the proton NV matrix element is very small and varies with all relevant 
parameters. 

We have also studied the effective nuclear interaction in the same CB used to calculate the 
weak interaction matrix elements. The dominant static part of the lowest order 2-body 
gives fairly accurate results for the spin, isospin and spin-isospin susceptibilities of nucleon 
matter. However, it is necessary to include at least 3-body effects to obtain the minimum 
in the Eo{p) of symmetric nuclear matter. The is much weaker than the bare v, and 
presummably can be used in perturbation theory formalism. 

All calculations of weak transition rates using effective interactions must in principle use 
the quenched matrix elements calculated in the same basis. We plan to calculate the weak 
interaction rates in nucleon matter using the present effective operators and interactions. 
To obtain more accurate predictions, it will be necessary to include > 3-body terms in the 
cluster expansion of the CB effective operators and interactions. 
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APPENDIX A: SECOND ORDER PERTURBATION THEORY 


Standard perturbation theory is applicable when the bare interaction Vij is weak. We 
then have H = Hq + Hj, Ho = T and 


Hi = J2v 

i<j 


ij ? 


(Al) 


Let |<hx) be the unperturbed FG state. The perturbed, normalized state up to second order 
is given by: 


|X) = ( 1 - - ^ |($y|iLr|$x)P 


2 


l*x> + ^ l*y> 

Yj^x 


(<hy|iL/|$x) 


^XY 


Ia> \{^y\Hi\^z) {^z\Hi\^x) las \{^y\Hi\^x) {^x\Hi\^x) \ ^ ^r,^ 

^ -- ^ ^ - 


E'^xz 


FP 

^XY 


FP 

^XY 


Y,Zj^X Y+X 

E^y = Tpg{X) — TpciY). In this approximation the Fermi matrix element is given by 


{E\Op\I), where <F/ and ^p are given by Eq. (|l^ . 

We are concerned only with two-body effects and therefore consider only the interactions 
Vaj in Hp The last two terms of the above |X) can be combined with the second by replacing 
the Vaj by an effective interaction; hence we will omit them. The direct terms of (FlOirlJ) 
can be written as: 

(oE Op{i)\I)direct = E.Q.d.a + E.l.d.j + E.2.d.j + E.2.d.a , (A3) 

i 

since E.O.d.j and E.l.d.a are zero. E.n.x.y is defined as in Section || with the exception that 
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n here refers to the order of Hj. We obtain: 


F.O.d.a 

F.l.d.j 


F.2.d.j 

F.2.d.a 


(kp|Oi7(a)|k„) = 1 

^(kp,h^|Oi.(j) 

hjv 

'^{'kp,hN\Vaj — 
hjv 

y^(kp,hjv|naj 


Q 


Eo-Ho 

Q 


hjv 

E 

hjv . 
1 


(kp, h^v \ Vaj 


Eq + uj — Ho 

Q 

Eo + UJ — Ho 

Q 


^aj I kn, h JY ^ 

- OF{j)\kn,I in) 


OfU) ^ „ r;aj|k„,h^) 
-C/n — -nr 


0 ~ -'-'0 


Eo + UJ — Hn 


0_p(a) 


Q 


Eo-Ho 


'^aj I knj hjv) 


-(kp,hjv|Oir(a) Vaj rr F ^ rr Wi|k„,h^) 

/ -no — -Ho -no — -Ho 


l/v n I ^ 


Q 


Eo + UJ - Ho Eo + uj - Hi 


Vaj OF(a)|kn,h7v) 


(A4) 


(AS) 

(A6) 


(A7) 


where Eo = e(k„) + e(h^r), uj = e(kp) — e(k„), Q is the projection operator to ensure Pauli 
exclusion in intermediate states, and Hat are any occupied proton or neutron states. We use 
e(k) to denote single particle energies; when Ho = T, e(k) = /2m. 

In order to make a connection with the correlated basis theory, we see that in perturbation 
theory the unnormalized two-body wave function is given by: 


|vI/) = 



(A8) 


Comparing it with the correlated wave function (Eq. 

Q 


we can identify: 


(F„ - 1) 


(Bo - -ffo 




(A9) 


when the interaction is weak. In reality, Vij is strong and Eq. 


is not useful. The 


correlation operator is determined variationally and its uj dependence is neglected assuming 
that the average value of Eo — Ho is much larger. 

It can be verihed that all of the F.n.d.y terms in Sect. II are obtained by replacing the: 


in Eqs. A4 to A7 by {Faj 


Q 

Eo-Ho 


Vaj 


and Vaj 


Q 

Eo + UJ — Ho 


1), since F^ = F. 
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APPENDIX B: THE C- AND F-COEFFICIENTS 


The C-parts required to calculate the effective weak vector operators in CB are obtained 
as follows: Let X, Y, Z be operators of type: 


X = . 

p=l ,6 

The C-part of the product of operators is then given by: 

c{xYZ) = Y, Y ’ 

^, 9 = 1,6 r,s=l ,6 


where = 1, and the are given in Ref. The results are listed below. 

= {rf + ^{rf + ^{tf , 
c'd' = cf = if - i)r + ^rr + , 

Cf = (r-l)2 + 3(r)2 + 6(f)^ 

^00 ^ ^ 32(/*)2 + 6(r - i)r, 

= (/-)2 _ 3(r-)2 + + 6/V , 

cf = cf = (r - i)r - ^rr + 12 /v*" + 3(r - fr + srr 


(Bl) 


(B2) 


(B3) 

(B4) 

(B5) 

(B6) 

(B7) 

(B8) 


The (Ta and A* parts of a product of cTq ■ crj, Saj , cr^ and CTj operators is obtained by 
repeated use of the Pauli identity: 


cr • Bcr •C = BC + icrBxC 


(B9) 


to reduce it to terms linear in cr^, cTj. Terms linear in <Tj go to zero on summing over j. The 
remaining terms linear in cr^ are expressed in terms of the operators cr^ and A^ to obtain 
the following equations. 


^00,(T 

d,a 


(BIO) 

p01,a 

^ d,a 

= {p-Y)f -rr -2fr, 

(Bll) 

plljCr 

d,a 


(B12) 

rOO.n 
^ d,a 

= 4^/* + 2(/*)^ 

(B13) 

p0l,A 

= 2/V*" + 2/V + 2/V‘" , 

(B14) 

pll,A 
^ d,a 

= Arf^ + 2{ff^ , 

(B15) 
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T7’00,cr _ 

^d,j - 

_ 7^01,(T _ 

- ^d,j - 

T7ll,cr _ 

^d,j - 

7700,A 

^d,j - 

plO,A _ p01,A _ 

^d,j - ^d,j - 

plliA _ 

^d,j - 


2{r-i)r + A{rf-A{pf , 

(r - i)r + rr + 2 rr" - 2 ff ^, 

2rr + 2{rf-2{f^)\ 
2ir-i)f-2rf + 2iff , 

ir - 1)/*" + /V" - /"/*" - /V + 2/v*", 

2jrjir _ 2f<^rjtT 2{f^f . 


(B16) 

(BIT) 

(B18) 

(B19) 

(B20) 

(B21) 

(B22) 


^eT = (/' -1)' + 2(r - i)r + iry - , (b23) 

F!!a^ = ir - rr - ir - rr + rr + rr - rr , (b24) 

Fir = rrr - 1 ) + nr - 1 ) + rr - rr - rr , (b25) 

Fir = rr + ir? + ir? - rr ?, (b 26 ) 

= 2 (r - rr + ‘^rr+rr ?, (b27) 

T"e,a’^ = rr - rr - ‘^rr+rr+rr, (b 28 ) 

Fi?^ = -2rr +4/*v"+2/v*+ rr ', (b29) 

Fi?^ = 2 rr+ 2 rr+rr ?, (b3o) 

Fir = ir - rr + ir - 1 ) - rr? + ir?, (b3i) 

Fir = ir - rr + rr - rr - rr - rr + rr, (b32) 

Fir = rir - 1 ) + nr - 1 ) - rr' + rr + rr , (b33) 

Fir = + ir? - rr? + ir ?, (b34) 

Fir = rr - 1 )+ rf+rr?, (b35) 

Fir = rr + rr - rr+rr , (b36) 

f^a ^ 2n{r -1)+ rr' - ^rr+rr , (b37) 

Fiy = 2 r^r+‘^rr+rr? ■ (b38) 
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I 


P{a,Po) 

P{a, Ipo) 

n t 

0.92 

0.89 

0.87 

n i 

0.02 

0.03 

0.03 

p T 

0.02 

0.03 

0.03 

p i 

0.04 

0.05 

0.07 


TABLE I: Correlated Basis probabilities for the active quasi-nucleon a to be | and A J, in the 
initial state for p = ^, 1, |/9o and Xp = 0.5. The listed values include contributions of 1- and 
2-body direct terms. 


P 

\{OF)?{q = Q) 

Otp 

PGTiq = 0) 

acT 

0.08 

0.80 

-0.094 

0.76 

0.259 

0.16 

0.81 

-0.075 

0.75 

0.060 

0.24 

0.86 

-0.083 

0.78 

0.041 


TABLE II: Quadratic fit to rjp and r]Qp at small q. 


P 

Cp{p) 

Cn{p) 

0.08 

1.39 

-1.29 

0.16 

1.53 

-1.46 

0.24 

1.40 

-1.38 


TABLE III: Linear fit to 

VNA for 

Xp < O.E 

t at small q. 

Density (fm“^) 

0.08 

0.16 

0.24 

1-b Tfg 

13.9 

22.1 

29.1 

2-b all 

-25.9 

-43.7 

-56.2 

3-b static 

4.9 

10.9 

19.1 

3-b LS + >4-b all 

-2.2 

-1.7 

0.8 

{Eq - Ey) 

-0.6 

-1.8 

-3.3 

Calculated Eq 

- 9.9 

- 14.2 

-10.6 

Empirical Eq 

- 12.1 

- 16.0 

-12.9 


TABLE IV: Contributions to the ground state energy of SNM from Argonne Vij and Urbana Vijk 
in MeV per nucleon 
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F.l.e.a.N F.l.e.a.D 



F.2.d.a.N F.2.d.a.D 



F.2.e.a.N F.2.e.a.D 

FIG. 1: Diagrams illustrating all of the one and two-body terms contributing to the Fermi CBME. 
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FIG. 4: Correlation dependence of rjp for k]\f = kpN and p = po- The dashed line shows results with 
juT _ jtT _ addition, /'^ = 1 for the dash-dot line. The dotted line shows | ^ F.n.x.ap 

when = 0. The solid line gives the full result. 
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FIG. 5; Dependence of rjp on the initial (fen) and final (kp) momenta for p = pQ. Each set contains 
six lines depicting the results for kn = (.5, .75, l)A:_Fn, and Xp = 0.3 and 0.5 for the indicated value 









q (fm' ) 


FIG. 8; Correlation dependence of rjcT for for = kpN and p = pQ. The solid line is for pcT 
with the full F. The dashed line shows results with = 0, and in addition, f^ = \ for the 

dash-dot line. The dotted line shows | a, GT.n.x.ap when = 0. 





FIG. 9: Dependence of the rjcT on the initial (kn) and final (kp) momenta for p = pQ. Each set 
contains six lines depicting the results for kn = (.5, .75, l)kFn-, and Xp = 0.3 and 0.5 for the 
indicated value of kr,- 



Proton NV CBME / FGME 









q (fm‘ ) 


FIG. 12; Correlation dependence of the neutron r/jvy for ki = kf = kpn and p = pQ. The dashed 
line shows results with = 0, and in addition, f^ = l for the dash-dot line. The solid line 

gives the full result. 






Proton 



FIG. 14; Proton r]i\fA as a function of q and proton fraction Xp for ki = kj = kpp- The solid, 
dashed, dotted and dash-dot lines show results for Xp = 0.5, 0.4, 0.3, and 0.2. The stars are results 
for TjGT at Xp = 0.5. 












FIG. 15: E(jt{p) (upper set), Er{p) (middle set) and Ea{p) (lower set) of symmetric nuclear matter. 
In each set, the upper most curves are results using Eij ioT p = ^po, the middle for p = po, and the 
lowest for p = |/9o- Solid lines show the results for and the dashed lines . Stars denote 
values obtained for Et{p) from variational calculations |27]. 
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FIG. 16: Ea{p) for pure neutron matter. The solid line shows results obtained using and the 
dashed for the . The results obtained with Fij for p = |, 1, |/9o are essentially indistinguish¬ 
able. Stars denote values obtained for (p) from quantum Monte Carlo calculations The 

dash-dot line is the Fermi-gas Efj(p). 
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FIG. 17: E(j{p) for symmetric nuclear matter (lower set of curves) and pure nuetron matter (upper 
set of curves). In each set, the upper most curves are results using Fij for p = ^po, the middle 
for p = pq, and the lowest for p = Solid lines show the results for and the dashed lines 
. Stars denote values obtained for Eo{p) from variational calculations |27]. 
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FIG. 18: The Fourier transform of the central and cr* • aj components of using Fij obtained 
at p = 1, |/9o are shown by dotted, solid, and dash-dot lines respectively. The dashed line shows 

the Fourier transform of the corresponding bare interaction. 
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q (fm‘' ) 


FIG. 19: The Fourier transform of the r* • Tj and cTj • (Tjn • Tj components of using Fij 

obtained at p = 1, |/9o are shown by dotted, solid, and dash-dot lines respectively. The dashed 

line shows the Fourier transform of the corresponding bare interaction. 
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FIG. 20; The Fourier transform of the Sij and r* • TjSij components of using Fij obtained 

at p = 1, |/9o are shown by dotted, solid, and dash-dot lines respectively. The dashed line shows 

the Fourier transform of the corresponding bare interaction. 
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